function [ ] = plotDirectionField( pdir, pev, domain2D )
%PLOTDIRECTIONFIELD Summary of this function goes here
%   Detailed explanation goes here

px = [domain2D.vx];%; rx];
py = [domain2D.vy];%; ry];

n = size(px, 1);

dir = zeros(n, 2);

pev(pev<39) = 0;

for i = 1 : n   
   dir(i, :) = pdir(i, :) * pev(i);    
end


figure;

quiver(px, py, dir(:,1), dir(:,2));

% axis ( [ -0.4, 0.4, -1, 1] );
% caxis( [ 0, 1.7] );
axis square;

end

